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Abstract 



We calculate the 0{as) virtual corrections to the matrix element for b —>■ sj, taking into account 
the contributions of the four-Fermi operator O2 and the electromagnetic and color dipole-tj^e 
operators. The results are combined with existing 0{as) Bremsstrahlung corrections in order to 
obtain the relevant inclusive rate. The new result drastically reduces the large scale dependence 
of the leading logarithmic approximation. It implies that a very accurate prediction for the 
branching ratio for B Xg'^ will become possible once also the corrections to the Wilson 
coefficients are available. 
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1 Introduction 



In the Standard model (SM) flavor-changing neutral currents only arise at the one-loop level. 
The corresponding B-meson decays are therefore particularly sensitive to 'new physics'; but 
also within the Standard model framework, they can be used to constrain several Cabibbo- 
Kobayashi-Maskawa matrix elements involving the top-quark. For both reasons precise exper- 
imental and theoretical work is of great importance. 

B — > K*'-f is the first rare B decay mode, which has been measured in 1993 by the CLEO 
collaboration and recently also the first measurement of the inclusive photon energy spectrum 
and the branching ratio in the decay B Xg +'~f was reported . In contrast to the exclusive 
channels, the inclusive mode allows a less model-dependent comparison with theory, because 
no specific bound state model is needed for the final state. 

This data is in a good agreement with the SM-based theoretical computations presented in 
P, given that large uncertainties exist in both experimental and theoretical results. In 
particular, the measured branching ratio BR{B — > X^'y) = (2.32 ± 0.67) x 10~^ overlaps 
with the SM-based estimates in [|, || and in @, 0] • 

As the experiments are becoming more precise in the near future, also the calculations 
must be refined in order to allow quantitative statements about new physics or standard model 
parameters. 

It is well known that QCD corrections to the decay rate for b —>■ s'j bring in large logarithms 
of the form a^{mw) ^og"^{mb/M), where M = rrit or rrtw and m < n (with n = 0,1,2,...). 
These large terms can be resummed by renormalization group techniques. At present, only the 
leading logarithmic corrections (i.e. m = n) have been calculated systematically. In this work 
we include one class of next-to-leading corrections which we describe in more detail below. 

The calculations are most easily done in the framework of an effective theory which is 
obtained by integrating out the top quark and the VT-boson in the standard model and other 
heavy particles in extensions thereof. A complete set of dimension-6 operators relevant for the 
process 6 — >■ 57 (and b — >■ s'jg) is contained in the effective Hamiltonian P] 

^e//(fe-37) = -^A,EQ(/^)0.(/i) , (1) 

V2 

where Gp is the Fermi constant coupling constant and Cj{fi) are the Wilson coefficients eval- 
uated at the scale /x, and Xt = VttV^l with Vij being the CKM matrix elements. The operators 
Oj read 

01 = {cLpl^bLa) {sLal,,CLp) , 

02 = {cLaYbLa) (sl/^T^CL/?) , 

O7 = (e/167r2) a'^'^ (m;,(/i)/? + m,(/i)L) 6„ F^,, 

Og = {gs|lQ^^^)-So.<y^''{m,{^^)R + m,{^J)L){X%/2)bpG^^, , 

(2) 

where e and Qs are the electromagnetic and the strong coupling constants, respectively. In the 
magnetic moment type operators 0^ and Og; Ftiv and Gj^^ are the electromagnetic and the 
gluonic field strength tensors, respectively and L = (1 — 75)/2 and R = {1 + 75)/2 stand for 
the left and right-handed projection operators. We note here that the explicit mass factors in 
O7 and Og are the running quark masses. We did not give explicitly the four-Fermi operators 
Os-Og in eq. (0), because they have small and negligible Wilson coefficients. 
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To leading logarithmic precision, it is consistent to perform the matching of the effective 
and full theory without taking into account QCD-corrections and to calculate the anomalous 
dimension matrix (8 x 8) to order [|10l- The corresponding leading logarithmic Wilson 
coefficients are given explicitly in |ll|]. The leading logarithmic contribution to the decay 
matrix element is then obtained by calculating the tree-level matrix element of the operator 
C7O7 and the one-loop matrix elements of the four-Fermi operators CjOj {i = 1, ...,6). In the 
NDR scheme which we will use in this paper, the latter are non-zero only for i = 5,6. Their 
effect can be absorbed into a redefinition of C7 C^-^-^ ^ 

C!^ff = C7 + QdC, + 3QdC, . (3) 

Since the first order calculations contain large scale uncertainties, it is important to take into 
account the next-to-leading order corrections. A complete next-to-leading calculation contains 
two classes of improvements: First, the Wilson coefficients are required to next-leading order 
at the scale /i ~ nib- This requires the matching with the full theory (at /i = mw) at the 0{as) 
level and the renormalization group equation has to be solved using the anomalous dimension 
matrix calculated up to order a^. Second, the real and virtual 0{as) corrections for the matrix 
element (at scale fi ~ rrib) must be evaluated. 

The higher order matching has been calculated in ref. [T^ and work on the Wilson 
coefficients is in progress. In the present paper we complete the second step. While the 
Bremsstrahlung corrections have been worked out ^ [13], Q in order to get a non-trivial 
photon energy spectrum at the partonic level for B — > Xg'y, the virtual corrections to 6 — S7 
have not been completely known so far. Only those virtual diagrams which are needed to can- 
cel the infrared sigularities generated by the Bremsstrahlung diagram were calculated. In the 
present paper we evaluate all the additional virtual correction, neglecting, as mentioned the 
contributions of the small operators O^-Oe. These new contributions substantially reduce the 
strong scale dependence of the leading calculation. 

In the following we thus consider the virtual 0{as) corrections to 6 ^ 57 due to the four- 
Fermi operator O2 and the magnetic operators O7 (which has already been calculated in the 
literature) and Og (which is new); note that the operator Oi does not contribute to the matrix 
elements for b ^ sj and b — > s'jg because of its color structure. As the corrections to O7 and 
Os are one-loop diagrams, they are relatively easy to work out. In contrast, the corrections to 
O2, involve two- loop diagrams. 



2 Virtual Corrections to O2, O7 and Os 

Since the virtual and Bremsstrahlung corrections to the matrix elements are only one (well- 
defined) part of the whole next-to-leading program, one expects that this contribution alone will 
depend on the renormalization scheme used. Even within the modified minimal subtraction 
scheme (MS) used here, two different "prescriptions" how to treat 75, will lead to different 
answers. Since previous calculations of the Bremsstrahlung diagrams have been done in the 
NDR scheme where also the leading logarithmic Wilson coefficients are available, we also use it 
here. A discussion of the results in the 't Hooft-Veltman scheme (HV) |]15| is presented in ref. 



16| . While the one-loop (a^) matrix element of the operator O2 vanishes, we must consider 



several two-loop contributions (aj) shown in Fig. |I]. The diagrams in Fig. |I]are grouped in such 



For the analogous b ^ sg transition, the effects of the four-Fermi operators can be absorbed by the shift 
C^^^ = Cs + C5. 

This omission will be a source of a slight scheme and scale dependence of the next-to- leading order result. 
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Figure 1: Non-vanishing two-loop diagrams associated with the operator O2. The fermions {b, 
s and c quark) are represented by sohd hues. The wavy (dashed) hne represents the photon 
(gluon) . 

a way that the sum of each hne is gauge invariant. The corresponding two-loop integrals are 
calculated by standard Feynman parameter technique. The heart of our procedure which we 
explain in detail in ref. |jl6|] is to represent the rather complicated denominators in the remaining 
Feynman parameter integrals as complex Mellin-Barnes integrals O, |18|, |19|, EH] . After inserting 



this representation and interchanging the order of integration, the Feynman parameter integrals 
are reduced to well-known Euler beta-functions. Finally, the residue theorem allows to represent 
the remaining complex integal as the sum over the residues taken at the pole positions of certain 
beta- and gamma-functions; this naturally leads to an expansion in the ratio z = {rric/mb)'^, 
which numerically is about z = 0.1. It is, however, not a Taylor series because it also involves 
logarithms of z. A generic diagram which we denote by D has the form 



D = co + 



n,m 



^^"l0g" 



mt. 



mt 



(4) 



where the coefficients Cq and c„m are independent of z. The power n in eq. (^) is in general a 
natural multiple of 1/2 and m is a natural number including 0. In the explicit calculation, the 
lowest n turns out to be n = 1. This implies the important fact that the limit rric —>■ exists; 
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thus, there cannot be large logarithms (if the charm quark is replaced by the up quark with its 
small mass) in these diagrams. As we show in |]16[, the power m of the logarithm is bounded 
by 4 independently of the value of n. We also proved that the expansion converges, at least for 
z < 1/4. In our results, we have retained all terms up to n = 3. The numerical result truncated 
at n = 3 differs from the result for n = 2 by only about 1%. 

The final result for the dimensionally regularized matrix element M2 of the operator O2 
which represents the sum of all two-loop diagrams in [I| is (The individual sets of diagrams are 
given in ref. 
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27 



X — CF{s'y\07\b)tree 
TT 



(5) 



Note that we introduced the renormalization scale in the form /i^exp(7£;)/(47r) which is conve- 
nient for MS subtraction. Here, z = {mc/mf,)'^ and L = log(2;) and e is the ultraviolet regulator 
in the convention, where c? = 4 — 2e. The symbol ( in eq. (^ denotes the Riemann Zeta func- 
tion, with C(3) ~ 1.2021; Cp = 4/3 is a color factor and the matrix element {s''y\OY\b)tree is the 
0{a^) tree level matrix element of the operator O7. In eq. (§) we have inserted the numerical 
values for the charges of the up and down type quarks {Qu = 2/3 and Qd = —1/3). 

There are also counterterms to be included. As we are interested in contributions to b ^ sj 
which are proportional to C2, we have to take, in addition to the two- loop matrix elements of 
C2O2, also the one-loop matrix elements of the four Fermi operators C2SZ2jOj (j = 1, 6) and 
the tree level contribution of the magnetic operator In the NDR scheme the only 

non- vanishing contributions to 6 — *■ 57 come from j = 5, 6, 7. The operator renormalization 
constants Zij can be extracted from the literature in the context of the leading order 
anomalous dimension matrix: 



SZ25 
SZ27 



a. 



487re 

a. 



bZ' 



26 



167re 



C, 



loTie 9 



Defining 

M2, = {s^\5Z2,0,\b) , 
we find the following contributions to the matrix elements 



(6) 
(7) 



-2e 



M25 
M26 



3as I (mb 

— — L^d<^F - — 
167r e V 



'2e 



{si\Or\b) 



tree 



4 



as ("iQuCp QdCF\ 1 , 1^ 



(8) 



We note that there is no one-loop contribution to 6 — > 57 from the counterterm proportional 
to (6*2^0^2), where the evanescent operator OI2 (see e.g. the last ref. in [1^) reads 



(7^ ^ l\lPlvlp\ 



Oo 



(9) 



Adding the two- loop expression M2 (eq.(^)) and the counterterms (eq. (|D), we find the renor- 
malized result, which can be written as 



with 



M2""" 



{si\0,\h) 



tree 



416 



I2 log — + r2 



(10) 



(11) 



243 

+ 

+ 
+ 



{-833 + UAti'^z^'^ 

1728 - 1807r2 - 1296C(3) + (1296 - 3247r2)L + 1081.2 + 36L' 
648 + 727r2 + (432 - 2167r2)L + 36L' 
"-54 - 847r2 + 1092L - 756L^ 



(12) 



167r r 



45 



37r2 + 9L + 9L2 



z + 



Here, 3ftr2 and S5r2 denote the real and the imaginary part of r2, respectively. 

The virtual corrections associated with the operator O7 have been taken into account by 
Ali and Greub, 



see e.g. 
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These corrections contain infrared singularities and for 
nis = also collinear singularities. At the level of the decay rate, these singularities cancel if 
one adds the Bremsstahlung correction based on the matrix element squared of the operator 
O7. Denoting by Tjj the sum of the lowest order contribution, the virtual 0{as) correction and 
those parts of the Bremsstrahlung correction just mentioned, this finite quantity reads in the 
limit rris = 

a, /16 4772 rrib' 

1 + — 1^^ — +41og — 

3 /i ^ 



77 



37r V 3 



(14) 



where the lowest order contribution Tj^ is given by 



i 77 



327r4 



Ole 



(15) 



For later convenience, it is useful to define a 'modified' matrix element for h 57, in such 



a way that its square reproduces the result in eq. 
reads 



This modified matrix element M^°'^ 



Ml 



mod 



(16) 
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with 

^7 = ^ , r, = ^(4-vr2) . (17) 

Finally, we consider the contributions to 6 ^ S7 generated by the operator Og, i- e., the 
matrix element 

Ms = {s^\Os\b) . (18) 
The corresponding Feynman diagrams are shown in Fig. |^. The sum of all 6 diagrams in Fig.^ 



Ob 

— a X , X 



Ob 
X 9 



Figure 2: Contributions of Os to b ^ sj. The cross (x) denotes the possible place where the 
photon is emitted. 



yields 



QdCp as 



12 

33 + 27r^ + 241og(mb//i) - 6in 

e 



{si\Or\b) 



tree 



(19) 



12 TT 

There is also a contribution from a counterterm, reading 

Ms{counter) = dZ^j {s-i\0'j\b)tree , (20) 
where the renormalization constant 

(T 1 

6Zsr = -CFQd- (21) 
vr e 

can again be taken from the literature [|1^. We thus arrive at the renormalized result Mg: 



Mr = {S7\0,\b),ree ^ [^8 log y + Tg 



with 



32 



27 



-33 + 2n^ - Gin 



(22) 



(23) 



3 Impact on the branching ratio 

To summarize, we have calculated the virtual corrections to 6 — > 57 coming from the operators 
O2, O7 and Og. The contributions from the other operators in eq. (|l|) are either small or 
vanish. As discussed above, some of the Bremsstrahlung corrections to the operator O7 have 
been transferred into the matrix element M^°'^ for 6 — 57 in order to simplify the following 
discussion. While we have neglected those virtual corrections of the operators O^-Oq which are 
given by the analogous diagrams as shown in Fig. |1|, we took into account the non-vanishing 
diagrams of O5 and Oq where the gluon connects the external quark lines and the photon is 
radiated from the charm quark; these corrections are automatically taken into account when 
using Cy-^-^ instead of Cj. Since the Wilson coefficients of the omitted operators are about fifty 
times smaller than that of the leading one and since we expect that their matrix elements can 
be enhanced at most by color factors, it seems reasonable to neglect them. 
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We can now easily write down the amplitude A{b 57) for 6 — > 57 by summing the various 
contributions derived in the previous section. We follow closely the treatment of Buras et al. 
0, where the general structure of the next-to- leading order result is discussed in detail. We 
write 

A(6^S7) = D {s^\OM\b),ree (24) 



with D 



D = C^/^if.) + ^ (cr^'^'m log ^ + Cr^'^rA , (25) 



and where the quantities ii and are given for i = 2,7, 8 in eqs. (|rT],|T^, (|I^) and 
respectively. The first term, Cj^^{fi), on the rhs of eq. (p5D has to be taken up to next-to- 



leading logarithmic precision in order to get the full next-to- leading logarithmic result, whereas 
it is sufficient to use the leading logarithmic values of the other Wilson coefficients in eq. 
(|25|) . As the next-to-leading coefficient Cj^^ is not known yet, we replace it in the numerical 
evaluation by its leading logarithmic value Cf''^'^^ . The notation {s'')\0^{ij)\h)tTee in eq. (|2^) 



indicates that the explicit mf, factor in the operator O7 is the running mass taken at the scale 

Since the relevant scale for a h quark decay is expected to be ~ m^, we expand the matrix 
elements of the operators around /i = up to order 0{a^). The result is 

AC A 

A{h ^ S7) = D {s^\0,{m,)\h),,ee (26) 

D = a;"i,) + ^ (Cf )^^^(/.)7i?^^^ log ^ + ^^^r.) , (27) 

where the quantities '^^^^^^^ = + 8 di-j are just the entries of the (effective) leading order 
anomalous dimension matrix [^]. As also pointed out in this reference, the explicit logarithms 
of the form Q;s(mf,) log(mb//i) in eq. ( p?]) are cancelled by the /^-dependence of Cf'*^^^^). 
I Therefore the scale dependence is significantly reduced by including the virtual corrections 
calculated in this paper. 



with D 



The decay width F''*'^* which follows from A{b 57) in eq. ( 126|) reads 



^virt _ ^ b,pole G\\^ae m 

327r4 



F\D\^ , (28) 



where in fact we discard the term of 0{al) in \D\ . The factor F in eq. (pSf) is 



^ ^ / mb(/i = mb) \ ^ ^ ^ _ 8 Q,(mb) 

V ^b^poie J 3 TT ■ 

To get the inclusive decay width for h — > s^y^g), also the Bremsstrahlung corrections (except 
the part we have already absorbed) must be added. The contribution of the operators O2 and 
O7 have been calculated before by Ali and Greub M, recently also the complete set has been 



worked out H, |T3|, IT^- In the present work we neglect the small contribution of the operators 



O3 - Oq in analogy to the virtual corrections, where only O2, O7 and Og were considered also. 



^As we neglect the virtual correction of Os-Oq, there is a small left-over dependence. 
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The branching ratio BR(6 s'~f{g)) is then obtained by deviding the decay width F = r"''*+ 
Ybrems ^ _^ s'y{g) by the theoretical expression for the semileptonic width Tgi (including the 
well-known 0{as) radiative corrections |^) and by multiplying with the measured semileptonic 
branching ratio BR,/ = (10.4 ± 0.4)% 

In Fig. ^we present the result for the branching ratio for b s'~f{g) based on the contribu- 
tions discussed above. A rather crucial parameter is the ratio mc/mb] it enters both, b —>■ sj 
mainly through the virtual corrections of O2 and the semileptonic decay width through phase 
space. To estimate this ratio we put rrify^poie = 4.8 ± 0.15 GeV for the b quark pole mass; from 
the mass difference — rric = 3.40 GeV, which is known quite precisely through the l/mq 



expansion p3| , p^ , one then gets mc/mf, = 0.29 ± 0.02. In the curves we have used the central 
values for nih^poie and mc/mh. For the CKM matrix elements we put Vcb = Vts and Vtb = 1. In 
Fig. ^ we have plotted the calculated branching ratio as a function of the top quark mass rrit. 
The horizontal dotted curves show the CLEO limits for the branching ratio BR(i? ^s'l) PI- 
To illustrate the scale dependence of the branching ratio, we varied the scale /i between (mf,/2) 



160 170 



Figure 3: Branching ratio for b — > 57(5') based on the formulae in this section. The different 
curves are explained in the text. 

and (2mb). We considered two 'scenarios' which differ by higher order terms. First, we put 



the scale /i = mi, in the explicit a, factor in eq. (pTf) and in the correction to the semileptonic 
decay width, as it was also done by Buras et al. 0. The resulting /i dependence is shown by 
dash-dotted lines. Second, we retain the scale /i in the explicit a, factors; this leads to the solid 
curves in Fig. |^. In both cases the upper curve corresponds to = mb/2 and the lower curve 
to /i = 2mfe. We mention that the /z-band is larger in the second scenario and it is therefore 
safer to use this band to illustrate the remaining scale uncertainties. 

In Fig. ^ we show for comparison the leading logarithmic result for the branching ratio 
for 6 — ^ 57, based on the tree- level matrix element of the operator O7 and using the tree- level 
formula for the semileptonic decay width. Varying the scale /i in the same range as above, 
leads to the dash-dotted curves in Fig. ^. We have also plotted the result as is was available 
before the inclusion of the virtual corrections of O2 and Og (but with Bremsstrahlung and 
virtual corrections to O7 included) . This is reproduced by putting £2 = ^"2 = = ''^8 = in our 
formulae. As noticed in the literature 0, |I4[, the /i dependence in this case (solid lines) is even 
larger than in the leading logarithmic result. 

From the results in Fig. ^ it was relatively easy to read off a reasonable prediction for the 
branching ratio within a large error which was essentially determined by the /i dependence. 
In the improved calculation (Fig. |^) the fi dependence is significantly reduced, because all 
the logarithms of the form as{mb)\og{mb/ fi) are cancelled as discussed above. In the present 
situation it is, however, premature to extract a prediction for the branching ratio from Fig. |^. 
This only will be possible when also Cj-^'^ is known up to next-to-leading logarithmic precision. 
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Figure 4: Branching ratio for b — > s'~f{g). The leading logarithmic result is shown by the dash- 
dotted curve; The solid line shows the situation before the virtual contributions of O2 and Og. 
See text. 

But this result will, essentially, shift the narrow bands of Fig. ^, without broadening them 
significantly. Thus, a very precise prediction will become possible and renewed experimental 
efforts will be required. It is, however, rewarding to see that the next-to-leading result will lead 
to a strongly improved determination of the standard model parameters or to better limits to 
new physics. 
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